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The possibiUty of adiabatically cooling a spin-1 polar Bose gas to a spin-nematic phase is theo- 
retically discussed. The relation between the order parameter of the final spin-nematic phase and 
the starting temperature of the spinor Bose gas is obtained both using the mean-field approach for 
the high temperature and spin-wave approach for the low temperature. We find that there exists 
a good possibility to reach the spin-nematic ordering starting with spinor antiferromagnetic Bose 
gases. 

PACS numbers: 37.10. Jk,03.75.-b,75.25.-|-z 



A trapped Bose-Einstein condensate (BEC) loaded in 
periodic counter-propagating laser beams to form an op- 
tical lattice [l| becomes a very important technique to ex- 
plore the general properties of condensed-matter physics 
because of the flexible tunability of the parameters, for 
example, changing the atom-atom interaction strength 
using Feshbach resonances 0, Q or reducing the hop- 
ping constant by increasing laser intensity. Therefore 
the optical lattice can imitate the ion-lattice potential 
one encounters in solid-state materials without the dis- 
advantage that a real compound has - the impurity. The 
pureness of an optical lattice enables the exploration of 
the strong-correlated systems which is the main issue of 
the condensed-matter physics. 

Recently, a theoretical suggestion about the possibility 
to adiabatically cool a weakly interacting spin- 1/2 Fermi 
gas to an antiferromagnetic phase of Hubbard model has 
been made 0], however, experimentally it is difficult to be 
fulfilled due to the extremely low temperature one has to 
reach to cool Fermi gases. On the other hand, Bose gases 
with high spins in a lattice offer different rich phases, such 
as dimerized states 0, @ and nematic 0, Q orders. The 
increasing interests in spinor optical lattices provide a 
chance to discuss open problems in the theory of quan- 
tum magnetism, which only zero-temperature physics has 
been so far discussed Q . Especially the maturity of the 
technique in trapping spinor Bose gases and the large 
entropy contained in the high spin systems benefit the 
opportunity to adiabatically cool a spinor Bose gas to an 
exotic quantum magnetic state. In this Letter we discuss 
the finite-temperature properties of spin-nematic states 
and the possibility of achieving the spin-nematic state 
from a spin-1 Bose gas. For simplicity, we shall limit 
ourselves to a sample with no net magnetization. 

We start with a spin-1 Bose gas with antiferromag- 
netic interaction, for instance "^^Na [l3|, which ground 
state is a polar state with zero spin expectation value 
[m . Now we adiabatically launch periodic potentials us- 
ing conter-propagating laser beams to form an isotropic 
rectangular optical lattice in 3D. In the case of only one 
atom per each potential well, Bose-Hubbard model can 
effectively represent the mechanism for atoms in an op- 



tical lattice [I4I. The hopping constant t decays expo- 
nentially with the square root of the laser intensity and 
the atom-atom interaction strength Us depends on the 
total spin 5* = 0, 2, where 5* = 1 is excluded due to the 
identity of bosons with one atom per each well. In three 
dimensions, the on-site interaction is effectively a linear 
function of the s-wave scattering length: Us oc as- With 
a suitable large atom- atom coupling U/t 1, the Bose- 
Hubbard model can be replaced by an effective Hamilto- 
nian H = ^ij^, where 

H,, = J(S, • Sj) + K{S, ■ S,)2 + J-K, (1) 

with J,K defined as J = -2t^ /U2, K = -2t^/3U2 - 
At^/SUo and symbolized as the next-neighbor sites. 
For ^^Na atoms, for instance, since 02 ~ 52aB and 
ao ~ 46as with Bohr radius as, we have U2 > Uq > 0, 
therefore K < J < 0. In this case, \U2 - Uo\/Uo < 1, 
hence |J| - \K\ and \J - K\ < \Jl\K\. The pa- 
rameters can be mapped into a unit circle by defining 
7 = tan-i K/J. For K < J < 0, -37r/4 < 7 < -7r/2. 
The two limits depend on the atom-atom coupling: if 
IC/2 - C/ol < C/o < as for ^^Na, 7 -> -Sn/A. On the 
contrary, if U2 ^ Uq > 0, j ^ —7t/2. Since the scatter- 
ing lengths are tunable using Feshbach resonances, one 
can explore all parameter regime experimentally. 

A numerical calculation by Troy et. al. Q suggested 
that the ground state should be a spin-nematic state for 
an isotropic 3D lattice if — 37r/4 < 7 < —tt/2. A nematic 
state comes from the spontaneously symmetry breaking 
of 0(3) rotation symmetry, though conserving the time- 
reversal symmetry. Therefore, the expectation value of 
any spins vanishes: {{S") = 0, a = a;, y, z) due to the con- 
served time-reversal symmetry, while the quadrupole or- 
ders defined as Q"'^ = l/2((S'"S'^-hS'^S'")-4/3(S"^) have 
finite expectation values due to the broken symmetry of 
the 0(3) rotation. With these two features in mind, one 
can define a new set of basis: \x) = :^(— 11) + | — 1)), — 
^(|1)-H| -l)),|z) = |0), where = -1,0,1) are the 
eigenstates of S^. Using the new basis, it can be easily 
verified that S°'\l3) — ieaf3j\j) with a/Sj = x,y,z. The 
spin expectation value of a pure state \'ip) = ^ da\oi) is 
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given by {tp\S°'\^p) = ieaP'jdjsd*, hence, if da are all real, 
the spin expectation values vanish. Therefore, for ne- 
matic states, there exists a vector d — (d^, dy, d^) that 
breaks 0(3) symmetry. Without lost of generality, one 
can rotate the system to make the vector d as z direction. 
In this way, the five possible components of quadrupole 
moments with the normality constraint {d'^ + dy + dl = 1) 
for a spin-1 system can be reduced to a single component: 
dl - 1/3. 

The same argument described above for pure states can 
be applied to mix states. A density matrix of a single 
site is defined as p = Z]a,/3=x,y.2 The spin 

expectation values are thus given by 



(2) 



The spin ordering is therefore associated with the imagi- 
nary part of p. Since p is hermitian, the imaginary part of 
a density matrix denotes the antisymmetry part of it. If 
the density matrix is real and hence symmetric, the spin 
expectation values disappear, and the system conserves 
the time-reversal symmetry. A real symmetric matrix 
can always be diagonalized, thus with a suitable choice 
of axes: 



p = p^^\x){x\+pyy\y){y\+p^^\z){z\. 



(3) 



If, as in the T = state, the broken 0(3) symmetry is 
in the z direction, one has p^^ = p^^ ^ p^^. The ex- 
pectation values of the spin square thus have the form: 
((S'^)2) = {{Sy)"^) ^ {{S')'^). Using this broken sym- 
metry, the only nonzero element of the five quadrupole 
order parameters: S^Sf^ + S'"S'^(a ^ /?), (S'^)^ - (Sy)^ 
and (S*^)^ — 2/3 is the expectation value in z direction: 



= (4) 



3^ 3' 
Here we have used the normality condition: Trp — 1. 
With the definition: q = p"^ — i, the nematic system 
only depends on the parameter g, the normal vectors d 
and — d defined as unit vectors in ztz direction. 

The two-body Hamiltonian |T]) can be rewriten as 

62 = 



H^3 = eo/^^^ + e2/^f , where eo = -4tVC/o 



—At^ IU2 and P^p is the projection operator which 
projects the pair i and j into a total hyperfine spin S 
state. Since P/J*^ projects the states to 15* = 0), the oper- 
ator should be rotational invariant. From this symmetry 
observation, the only possible rotational invariant state 
is y z which has to be proportional to 

IS* = 0). Therefore 



(5) 



a..l3—x^,z 



(2) 

As for Pi;', one needs five orthonormal states which are 
also orthogonal to jS" = 0). Since the states correspond- 
ing to S* = 2 are symmetric under the coordinate ex- 
change i <-> J, the five orthonormal states can be written 



FIG. 1: (Color online) Order parameter of the nematic state 
vs. temperature in unit of ksTpf = 6\K\. Solid line is the 
mean-field result, while the red dashed line represents the 
spin- wave approach. The two approximations using linear 
and quadratic dispersion are shown in blue dotted line and 
green dotted-dashed line, respectively. 



down by constructing the second order polynomials cor- 
responding to a symmetric rank-2 tensor. One of the 
possible construction is \Iai3) = j + 



for (3, \Io) 



i\z)i\z)j- 

\y)^\y)j)■ 



\x)i\x)3~h\y)i\y)j) and 



Therefore 



Pip = T.\^){i\- 
{1} 



(6) 



The principle of the mean-field approximation is that 
one reduces a many-body problem to a one-body problem 
by singling out one site and averaging the other sites 
served as a mean field. A mean-field treatment for a spin- 
1 Bose system has also been done by Chen and Levy[l^. 
The effective Hamiltonian is thus a single site operator 
given by 



zTt,[p,H,,] ^ zTt, p.ieoPp^ +e2P,P) , (7) 



with the coordinate number z. For a rectangular three- 
dimensional lattice, z = 6. Using Eqs.Q, ^ and ([6]) 
one can obtain the effective Hamiltonian as 



eff 



a—x,y,z 



(8) 



The self-consistent equation for the only nonzero 



quadrupole order q 



which comes from the re- 



lation p^- 



= p-H^Ss/kBT 



/Tre-"""/'"''^ reads 



3q 
2 



exp(|i§f) + 2- 



(9) 



q = is always a solution, actually it represents the 
disordered phase = p**^ = p^^ = 1/3, correspond- 
ingly, the number of particle per site in the state 



3 



FIG. 2: (Color online) Order parameter of the nematic state 
vs. starting temperature of spinor Bose gases in unit of Tjjgg 
using the mean-field approach (solid line) and the spin-wave 
theory (red dashed line). For the red dashed line, the inter- 
acting spectrum for the initial Bose gas is used, while the blue 
dashed-dot line shows the result with an ideal Bose gas as an 
initial state. 



jm^ — —1,0,1) obeys no — n±i. This state conserves 
all symmetries and gives the equal spin expectation val- 
ues: ((5")2) = {{Sy)^) = ((5^)2) = (^2^/3 = 2/3 as the 
result indicated by Eq. (U). 

The solid line in Fig. [T] shows the mean-field result for 
the order parameter p^^ as a function of T/T]\j where 
T/v = 6|iir|. In this figure we found a first order phase 
transition with the transition temperature Tc/T^ ~ 0.36. 
Below Tc, the system has the broken 0(3) symmetry in 
the z direction, therefore p^^ > 1/3. Above Tc, p^^ = 1/3 
for all temperaure. The order parameter jumps discon- 
tinuously from the ordered state to the disordered state 
at the transition temperature Tc- We have to mention 
that the transition temperature predicted by the mean- 
field theory is not correct as the case for the Ising model 
14| . and the order parameter in the low temperature 
decreases exponentially, as shown in Fig [U which is an 
artifact of the mean-field approach. However, the mean- 
field solution gives us a qualitative description about the 
ordered-disordered phase transition. 

The entropy per unit site given by the definition s = 
— fcsTrplnp reads: 



= -kf 



lnp^^ + (l-p^^)ln 



1 



(10) 



The critical order parameter p^^ = 0.67 at Tc, which gives 
the critical entropy: Sc = 0.863 according to the Eq. pO)) . 
The value of the critical entropy is quite large due to the 
large entropy capacity for spin-1 systems. Above Tc, the 
entropy jumps to a constant: Sd = ks In 3. The number 3 
indicates the mean-field solution takes equal probability 
for each conserved 0(3) symmetry for all temperatures 
T>Tc. 

At low temperature, the mean-field theory fails be- 
cause it ignores long- wavelength fluctuations. This prob- 



lem has been solved using spin-wave theory by Chubukov 
by approximating the spin-up state |1) and the spin- 
down state | — 1) as two bosonic states and |0) as a vacuum 
state. This approximation works very well if the excita- 
tion number per each site is less than one. The resulting 
energy spectrum consists two degenerate branches 

efe = (-zJ)(l + <5)i/2 {(1 - -uk + 6{l + , 

(11) 

where 5 — tan 7 — 1 and the tight binding function 
Vk = iZ-^X^r e*'^ '"'' with the all next-neighbor vectors 
r^. In the long- wavelength limit and for (5 <C 1, the en- 
ergy spectrum is linear: Sk = ck where c — —zJy^2S/za 
with the lattice spacing a. Therefore these two energy 
branches correspond to two Goldstone modes due to the 
two remaining symmetries for a broken 0(3) symmetry. 

The order parameter corresponds to the density 
still remaining in the vacuum state (remember \z) — |0)). 
In the momentum space, p^^ can be separated into two 
terms: zero temperature depletion r]{d) and thermal de- 
pletion rjT{5,T): p^^- = 1 — 770 (5) — riT{S,T), where 



VoiS) 




S+il-iyk) 



V(i-i^fc)[25+(i-z.fc): 

S + jl- i^k) 



- 1 



1) 



(12) 



The red dashed line in Fig [T] shows p^^ obtained by us- 
ing the spin- wave approach with the parameters K/ J = 
1.1 ((5 = 0.1), which are very close to the parameters 
K/J = 1.089 for 23iVa . Contrary to the result ob- 
tained from the mean field, for which p^^ = 1 at zero 
temperature representing the classical nematic ground 
state, the quantum fluctuation of the spin wave gives 
770 — 0.009. In Figd] two approximated curves are also 
shown. In the regime T/Tj^ <C S, where the energy 
dispersion is linear, r]T = \/3(l + S)'^/{2VS){T/Tn)^, 
as shown with the blue dotted line. Beyond the lin- 
ear regime, for S < T/Tn < 1, £fe « \J\k'^a'^. In 
this case the zero temperature fluctuation is zero and 
r?T = 2[3(1 + (5)/27r]3/2((3/2)(r/TAr)3/2 as shown with 
green dotted-dashed line. We can see that the order pa- 
rameter of spin wave theory coincides with the result us- 
ing linear dispersion at the temperature T/Tn < 0.1 and 
then goes parallel to the approximation with a quadratic 
dispersion. At higher temperature the spin-wave the- 
ory overestimates the depletion to the excited states and 
hence the order parameter decreases faster than the mean 
field approach. This is because at higher temperature, 
more than one boson can occupy one site. This feature 
can not explain the spin statistics more, therefore the 
spin-wave theory fails for the high temperature regime. 

The initial spinor Bose gas has three modes of 
excitations [lT|: one corresponds to density excitations 
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with the energy dispersion eo = [(A:^/2m)^ + ncofc^]^/^, 
and the other two correspond to two spin modes: e2 = 
[(A;^/2m)^ +nc2fc^]^/^ with the mass of atoms m and the 
density of condensate n. The interaction strengths Cq 
and C2 is defined as follows: cq = 47r(ao + 2a2)/3m and 
C2 — 47r(a2 — ao) /3m. If one starts with a Bose gas at the 
temperature ksT ^ ucq for which the interaction is not 
important, the system is approximated as an ideal Bose 
gas. The adiabatic cooling ensures that the beginning 
state (Bose gas) and the final state (nematics) should 
have the same entropy 



T 



5C(5/2) 

'2C(3/2) V^bec 



3/2 



(13) 



where T^ec = 27rn2/3/[3C(3/2)]2/3m, ({x) = V"" 
is the Riemann zeta function and s is given by Eq. ljlOp . 
The solid line in Fig. [2] shows the mean-field order pa- 
rameter of the cooled nematic state versus the starting 
temperature T/T^^^ of the Bose gas using the mean-field 
approach. There exits a maximum temperature for the 
starting state such that 7^/7bec < 0.769. On the other 
hand, for Sc < s < Sd, 0.679 < T/T^^^ < 0.9. In this 
initial temperature regime, the system is cooled to a co- 
existence of the ordered and disordered phase with phase 
separation due to the first order phase transition. For 
T/T^gj, > 0.9, the nematic phase can therefore not be 
reached. 

For low temperature, the entropy of bosonic systems 
has the form 



k.<T 



- ln(e'=B^ - 1) 



kBT{e' 



1) 



(14) 



where a represents different Goldstone modes and E'^ is 
the corresponding excitation energy with momentum k. 
For the nematic state, E'^ = Sk for two degenerate modes, 
while for the initial spinor Bose gases E^ — eq for one 
density mode and E^ = €2 for two spin modes. By map- 
ping the initial BEG temperature to the cooled nematic 
temperature with equal entropy, we found the order pa- 
rameter as a function of the initial BEG temperature as 
the red dashed line shown in Fig [21 At low temperature, 
the order parameter decays much sooner than that for 
the mean field, and the zero temperature fluctuation ex- 
ists. For comparison we also show the order parameter 
initially with an ideal Bose gas as the blue dotted-dashed 
line in Fig. [S] There almost exists no difference for p^^ 
between the interacting spinor Bose gas and the ideal 
Bose gas as an initial state because the chemical poten- 
tial nco for the interacting gas is small and hence is easily 
smeared out by the temperature. 

Moreover, we consider a BEG trapped in a harmonic 
trap with a trap frequency cu due to the fact that in reality 
a BEG of ultracold atoms only exists in a magnetic or 
optical trap. The nematic exists with uniform K and J 



in the case that the chemical potential at the center of the 
trap is less than the atom-atom interaction strength Us 
and the entropy at the edges is not dominant. The second 
condition can be possibly fulfilled because the nematic 
state carries a lot of entropy. The starting entropy for 
the trapped gas has the form: 



4C(4) 

C(3) \Tr 



T 



bee 



(15) 



where T^Q^. — n^^^Lij/[3({3)]^^'^. The starting tempera- 
ture regime for cooling into a nematic state is T/T^^^ < 
0.62 if s < Sc is used. Therefore the final result is similar 
to Fig. [2] with a squeezed temperature regime. 

p"" can be directly measured since it is the number 
of particles per site in the d • S = state where the 
quantization axis is chosen along the a-direction. The 
ordered state is characterized by p^^ > 1/3 > p^'^'^yy\ 
For quantization along z-direction, the ordered state can 
also be identified by rig > 1/3 > n±i. On the contrary, 
the disordered state has isotropic densities. However, 
in reality, there also always exists a small magnetic field 
which breaks the rotation symmetry due to the quadratic 
Zeeman effect Therefore for T > Tc, is slightly 
larger than 1/3, so it is not entirely disordered. However, 
for a small enough quadratic Zeeman field the first order 
phase transition still exists. 

In summary, we have shown, by means of both mean- 
field and spin-wave approaches, thermodynamically the 
order parameter of the nematic state for a spin-1 Bose 
gas in an optical lattice. With the mean-field approach, 
we found the first-order transition temperature. The en- 
tropy containing in a nematic state is very large as a 
result of higher spin systems. Using the spin- wave ap- 
proach, we found the zero-temperature fiuctuation and 
power-law decay of the order parameter. We further- 
more discussed the possibility of cooling a spinor polar 
Bose gas to a spin-nematic magnetic state. We found 
that there always exists a large temperature range under 
which the spin-nematic state can be reached. 
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